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Abstract 

We study by Wilson numerical renormalization group the spectral properties of a two-orbital 
Anderson impurity model in the presence of an exchange splitting which follows either regular 
or inverted Hund's rules. The phase diagram contains a non-Fermi liquid fixed point separating 
a screened phase, where conventional Kondo effect occurs, from an unscreened one, where the 
exchange-splitting takes care of quenching the impurity degrees of freedom. On the Kondo screened 
side close to this fixed point the impurity density of states shows a narrow Kondo-peak on top of a 
broader resonance. This narrow peak transforms in the unscreened phase into a narrow pseudo-gap 
inside the broad resonance. Right at the fixed point only the latter survives. The fixed point is 
therefore identified by a jump of the density of states at the chemical potential. We also consider the 
effect of several particle-hole symmetry-breaking terms. We show that particle-hole perturbations 
which simply shift the orbital energies do not wash out the fixed point, unlike those perturbations 
which hybridize the two orbitals. Consequently the density-of-state jump at the chemical potential 
remains finite even away from particle-hole symmetry. In other words, the pseudo-gap stays pinned 
at the chemical potential, although it is partially filled in. We also discuss the relevance of these 
results for lattice models which map onto this Anderson impurity model in the limit of large 
lattice-coordination. Upon approaching the Mott metal-insulator transition, these lattice models 
necessarily enter a region with a local criticality which reflects the impurity non-Fermi liquid fixed 
point. However, unlike the impurity, the lattice can get rid of the single-impurity fixed-point 
instability by spontaneously developing bulk-coherent symmetry-broken phases, which we identify 
for different lattice models. 

PACS numbers: 71.30.-|-h, 71.10.-w, 72.15.Qm, 72.15.Ru 
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I. INTRODUCTION 



Non-Fermi liquid behavior may emerge in Anderson and Kondo impurity models for two 
distinct reasons. The first one is that, by construction, the conduction electrons may not be 
able to perfectly Kondo-screen the impurity degrees of freedom. This is realized for instance 
in multi-channel Kondo models^. 

The alternative route towards non-Fermi liquid behavior is the presence of an intra- 
impurity mechanism which splits the impurity degeneracy favoring a non-degenerate con- 
figuration. The Kondo exchange takes advantage of letting the impurity tunnel among all 
available electronic configurations. This quantum tunneling is hampered by any term which 
splits the degeneracy and tends to trap the impurity into a given state. Therefore either the 
Kondo exchange overwhelms the intra-impurity splitting mechanism, or vice-versa, which 
leads respectively to a Kondo-screened phase or an unscreened one. When none of the two 
effects prevails, a non trivial behavior may appear. This is actually what happens in the 
two S=l/2-impurity Kondo model in the presence of an antiferromagnetic direct exchange 
between the impurity-spina^*^. There it is known that, under particular circumstances-', an 
unstable non-Fermi liquid fixed point separates the Kondo-screened and unscreened regimes. 
Since this fixed point requires fine tuning of the model parameters, it is tempting to con- 
clude that it is of little physical relevance. In reality a similar competition may be at the 
heart of strongly-correlated electron lattice-models. Here the kinetic energy profits by the 
electrons hopping coherently through the whole lattice. On the contrary, the strong cor- 
relation tries to optimize on-site (atomic) energetics, thus opposing against the hopping. 
This may involve two energy scales. The higher one is the so-called Hubbard U, which 
tends to suppress on-site valence fluctuations. The lower one, let us call it J, governs the 
splitting among on-site electronic configurations at fixed charge. It may be controlled by 
the exchange-splitting, the crystal field, by local distortion modes or even by short-range 
inter-site correlations. When the lattice model is driven towards a Mott metal-insulator 
transition (MIT), either by increasing U or by doping at large U, it necessarily encounters a 
regime in which the coherent quasiparticle bandwidth Wqp is of the same order as J, which 
we expect is essentially unaffected by U as it just determines the multiplet splitting at fixed 
charge. Since coherent hopping tends to occupy more or less democratically all multiplets, 
it opposes against J. Out of this competition interesting physical properties may emerge. 
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just like in the Anderson impurity models we discussed before. The analogy between the 
impurity and the lattice models can be even put on firm grounds in the limit of large coordi- 
nation lattices through Dynamical Mean Field Theory (DMFT)5. In that limit it is possible 
to map the lattice model into an effective Anderson impurity model (AIM) subject to a 
self-consistency condition which relates the impurity Green's function to the hybridization 
with the conduction bath. The quasiparticle bandwidth of the lattice model transforms into 
the Kondo temperature of the AIM. Since approaching the MIT Wgp 0, the effective 
AIM is necessarily driven into the regime ~ J, where the competition among the two 
screening mechanisms may result into anomalous physical properties. Exactly this compe- 
tition was invoked by Ref. to explain the appearance of a superconducting pocket, later 
shown to have a hugely enhanced superconducting gap '', just before the MIT in a model for 
alkali doped fuUerenes. 

More recently we have demonstrated by Wilson numerical renormalization group and by 
bosonization that a two-fold orbitally degenerate AIM in the presence of inverted Hund's 
rules possesses a non-Fermi liquid unstable fixed point similar to the two-impurity Kondo 
model one^i. Because of the aforementioned reasons, any lattice model which maps by DMFT 
into the same AIM should necessarily meet this fixed point on the route towards a MIT. We 
argued that, unlike the single- impurity, those lattice models might spontaneously generate 
by the DMFT self-consistency conditions a bulk order parameter to get rid of the single- 
impurity fixed-point instability. Since the fixed point is unstable in different, particle- hole 
and particle-particle, channels, there exist in principle several competing bulk instabilities. 
We speculated that, in the absence of nesting or band-structure singularities, the most likely 
instability is towards superconductivity. These predictions have been just recently confirmed 
on a lattice model by a DMFT calculation^. In this paper we pursue the analysis of that 
AIM by uncovering the spectral behavior across the non- Fermi liquid fixed point. This is 
not only interesting for the AIM itself, being one of the few cases where non-Fermi liquid 
dynamical properties may be accessed, but also in the context of the DMFT mapping. The 
model is also sufficiently simple to allow for an analytical description for the spectral function 
which reproduces well the numerical results and provide new physical insights. Actually our 
model spectral function has been quite useful in guiding the analysis of the DMFT solution 
presented in Ref. 0. 

The paper is organized as follows. In Section |n] we describe the two-orbital AIM model. 
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In Section ITTTl we introduce three lattice models which map by DMFT onto the two-orbital 
AIM: (a) a two-band Hubbard model e ® E Jahn- Teller coupled to local phonons; (b) a 
two-band Hubbard model in the presence of single-ion anisotropy; (c) two-coupled Hubbard 
planes. In Section Hyl we review in more details the Wilson numerical renormalization group 
calculations of Ref. ^ and present a new analysis based on Fermi liquid theory, which we 
develop in the Appendix. The new results concerning the dynamical properties are presented 
in Section In Section IVII we extract from the numerical data an analytical expression 
of the impurity spectral function. The role of symmetry breaking terms in particle-hole 
channels is investigated in Section IVIII Conclusions are presented in Section IVIIII 



II. THE MODEL HAMILTONIAN 

The AIM Hamiltonian we consider is 

H = Hu + Hj + H, + Hr,y, (1) 

= I (n, - 2 + z/)^ + 2 J [{T^f + iryf] + 

kacj 

Here c|^^^ creates a conduction electron in the band a = 1,2 with momentum k, spin a and 
energy ek, measured with respect to the chemical potential. is the creation operator 
of an electron with spin a in the impurity orbital a = 1,2, while Ud = J2aa ^aa^aa is the 
impurity occupation number. We have defined the orbital pseudo-spin operators 

T' = lT.J2^aa<Ao.^ (2) 

a a=l,2 

where i = x,y,z and r*'s are the Pauli matrices in the orbital space. We further assume that 
the conduction band density of states is symmetric with respect to the chemical potential, 
set equal to zero, so that the behavior of the Hamiltonian under a particle-hole symmetry 
transformation is controlled by the parameter in (^. For the time being we will take u = 0, 
which implies that the Hamiltonian is particle-hole symmetric. Afterwards we will release 
this constraint. The model without the impurity exchange coupling J is SU(4) invariant. 
A finite J lowers the SU(4) symmetry down to SU(2)spj„x 0{2)orbit- In this case the total 
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charge, the total spin and the total z-component of the pseudospin are the only conserved 
quantities. 

It is convenient to start our analysis by the spectrum of the isolated impurity, Vd = 0. 
The impurity eigenstates, \n, S, S^,T,T^), can be labeled by the occupation number n, the 
spin S, pseudospin T and their z-components, and T^, respectively, with energies 

E{n, S, S\ T, T') = ^{n-2f + 2 J [T(T + 1) - [T'f] . (3) 

We assume U ^ |J|, so that the impurity ground state with u = has n = 2. In this 
case the only configurations allowed by Pauli principle are a spin-triplet pseudo-spin-singlet, 
S = 1 and T = 0, 

|2,1, +1,0,0) = 4t4t 

12,1,0,0,0) = j.(^d\^dl-dl^dl^ |0), (4) 

12,1,-1,0,0) = 4^4J0), 
and a spin-singlet pseudo-spin-triplet, S = and T = 1. The latter is split by J into a 
singlet with = 0, 

|2, 0, 0, 1,0) = -^ (4^4 + 44^) |0), (5) 
and a doublet with = ±1, 

|2,0,0,1,+1) = JO), 
|2, 0,0, 1,-1) = dl^dl^ |0). 

If J > 0, the lowest energy configuration is the spin-triplet, S = 1 and T = 0, which 
corresponds to the conventional Hund's rules. On the contrary, for J < 0, the isolated 
impurity ground state is the singlet (0) with quantum numbers 5* = 0, T = 1 and = 0. 
We postpone to the following Section a discussion about physical realization of such inverted 
Hund's rules. 

A finite hybridization, Vd 7^ 0, induces valence fluctuations within the impurity, which 
are controlled by the energy scale (hybridization width) 

Ao = vr V,' p„ (7) 

with pc the conduction electron density of states (DOS) at the chemical potential per spin 
and band. These fluctuations are suppressed by a strong repulsion [/ ^ Aq, which we 
assume throughout this work. Although all our calculations refer to the AIM (0), it is more 
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insightful to discuss some physical properties in terms of the effective Kondo model which 
describes the low-energy behavior when U ^ Aq: 



^eff = Hj + Hc + Hk, (8) 
where Hj and He have been defined in (0) and the Kondo exchange 



Hk — Jk 



S-S + f-f + A WijW, 



i,j=x,y,z 



(9) 



with 



Here S, defined by 



Jk = 2V^/U. (10) 



a a/3 



T, which we introduced in Eq. and Wi 
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ah a/3 



are impurity spin, pseudo-spin and spin-orbital operators, respectively, while 5, T and 
Wjj are the corresponding conduction electron density operators at the impurity site. The 
impurity operators in Q act only in the subspace with two electrons occupying the d- 
orbitals, which, as we showed, includes six states. The Kondo model (jH)) contains two 
competing mechanisms which tend to freeze the left-over impurity degrees of freedom: (i) 
the Kondo exchange, with its associated energy scale, the Kondo temperature Tk\ (ii) the 
intra-impurity exchange splitting J. As we already mentioned, the Kondo exchange © gains 
energy by letting the impurity tunnel coherently among all available six configurations, but 
it is hampered by J which instead tends to trap the impurity into a well defined state. 

If J ^ Tk > 0, the positive exchange splitting dominates and the impurity is essentially 
frozen into the lowest energy spin-triplet configuration. The Kondo exchange projected onto 
the triplet sub-space (jH) is simply Hk = Jk S ■ <S, describing a standard S = 1 two-channel 
Kondo effect. This is known to be perfectly screened at low energyi^i^i, yielding a scattering 
phase shift 5 = 7r/2 in each spin and orbital channel. 

On the contrary, if J ^ —Tk < 0, the impurity gets trapped into the S = 0, T = 1 
and = configuration, Eq. (jSJ. Since (0) is non degenerate, the Kondo exchange is 



un-effective, so that asymptotically the impurity decouples from the conduction bath. This 
implies a low energy phase shift 6 = 0. The main question which we try to adress is how 
the model moves across the two limiting cases. 

As it was pointed out in Ref. 0, this behavior is parallel to the two S = 1/2 impurity 
Kondo model (2IKM) in the presence of a direct exchange between the impurity spins^iM. 
In that case, if the two spins are strongly ferromagnetically coupled, the model reduces to an 
5=1 two-channel Kondo model, while, if they are strongly antiferromagnetically coupled, 
the two spins bind together into a singlet and decouple from the conduction electrons, 
exactly as in our model. The two channels correspond in the 2IKM to the symmetric and 
antisymmetric combinations of the even and odd scattering channels with respect to the 
midpoint between the impurities. It was demonstrated by Ref. that, provided a peculiar 
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e-hole symmetry holds, the non-Fermi liquid unstable fixed point (UFP) found by 
separates the Kondo screened and unscreened regimes. In particular it was shown 



that while a particle-hole symmetry breaking term 

aa k,aa 

does not wash out the UFP, the latter is instead destabilized by the perturbation 

a 

-E^kcLck2a + ^-C-- (12) 

Translated into our two-orbital language, the dangerous symmetry which needs to be pre- 
served is just the 0(2)or{,j( orbital symmetry. Therefore, unlike in the 2IKM, where the 
two scattering channels are generically not degenerate, in our case the instability towards 
0{2)orbit symmetry breaking does correspond to a physical instability. Hence, if orbital sym- 
metry is unbroken, we do expect to find an UFP in our model, with similar properties as in 
the 2IKM. We notice that, in spite of the analogies, our model has a larger impurity Hilbert 
space than the 2IKM. In fact the S" = 0, T = 1 and = ±1 doublet of Eq. © is absent 
in the 2IKM, where it would correspond to doubly occupied impurities (the labels 1 and 2 
for the (i-orbitals translate in the 2IKM into the two one-orbital impurities). Yet we can 
perturb our Hamiltonian by adding to of the term 

Hg = G{T^)\ (13) 



with G > 0, which raises the energy of the doublet. If G ^ T^-, the doublet effectively 
decouples from the low energy sector, and our model should become equivalent to the 2IKM. 
In Section IIVI we show that indeed by increasing G our UFP smoothly transforms into the 
2IKM one. 

III. PHYSICAL REALIZATIONS 

As we emphasized in the Introduction, our interest in model (0) plus eventually 
is mainly motivated by its possible relevance for lattice models. In reality a formal cor- 
respondence bewteen single-impurity and lattice models holds strictly only in the limit of 
large lattice-coordination. Nevertheless we believe that this correspondence, at least close 
to a Mott transition, may remain valid even beyond that limit, making the single- impurity 
analysis of much broader interest. Therefore, although inversion of Hund's rules may indeed 
occur in realistic AIM's or in artificially designed quantum dot devices, here we rather focus 
on lattice models which map within DMFT into our AIM. 

A. Two-band Hubbard model in the presence of an e (8) -E Jahn- Teller coupling 

Let us start by considering a two-band Hubbard model in which each site is Jahn- Teller 
coupled to a doubly degenerate phonon. The Hamiltonian reads 

* = -4jt EE + H 

* a=l (7 <ij> 

i i 

+ Y E E (^'^ + pI) - 3 E + g., rn . (14) 

i a=x,y i 

Here —t/ ^/z is the hopping matrix element bewteen one site and its 2;-neighbors and Jh > 
is a conventional Hund's exchange, and qiy are the phonon coordinates at site i, and 
Piy their conjugate momenta, ujq the phonon frequency and g the Jahn- Teller coupling. The 
latter gives rise to a retarded electron-electron interaction whose Fourier transform is 
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If the phonon frequency ujq is much larger than the quasiparticle bandwidth we can safely 
neglect the dependence at low energy, so that the phonon-mediated interaction becomes 
un-retarded and given by 

LUq '—^ L J 

I 

Within DMFT the Hamiltonian maps in the large z-limit onto the same AIM model as in 
Eq. dH) with 

J=Ju~^. (15) 

which may be either positive or negative. The case with J < as well as the starting model 
realistically including phonons have been recently studied by DMPT^iil. 

B. Two-band Hubbard model with single-ion anisotropy 

Another realization which may also be physically relevant is the following lattice model: 

^ a=l <j <ij> 

i i 

+DY,{Stf. (16) 

i 

For Jh > and D ^ this model describes a two-band Hubbard model with conventional 
Hund's rules, favoring a spin-triplet two-electron configuration, in the presence of a single- 
ion anisotropy which splits the spin-triplet into a singlet with = and a doublet with 
= ±1. If D > 0, the 5*^ = configuration is favored. Upon interchanging S ^ T, this 
model maps in the z ^ oo limit onto with 

' = (17) 

G = D- 2Jh. 
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C. Two coupled Hubbard planes 



Finally let us consider two coupled single-band Hubbard planes described by the Hamil- 
tonian 

^ = -4j t E E + + 1 E - 

* a=l (7 <ij> a,i 

+ ^ JSi,-4 + v^K-i)K-i), (18) 

i 

where a = 1,2 labels the two planes and —tj^ is the in-plane hopping between one site 
and its 2;- neighbors. In the limit z ^ 00^ this model maps by DMFT onto an AIM self- 
consistently coupled to a batb^. The relations between the interaction parameters of the 
AIM, (P) plus (Uni), and those of ^ are given in Table HI 

TABLE I: Mapping between the AIM interaction parameters and the two Hubbard plane ones. 



AIM 


Two Hubbard planes 


U 


i(C7 + F)-ij 


J 




G 


U-V + \J 



In reality it is more interesting to consider the model ()18j) with J = V = Q but in the 
presence of an inter-plane hopping 

-t± (^^-^2.. + H.c) . (19) 
lo- 
in the limit of large lattice-coordination, this model maps close to the MIT onto a two- 
orbital AIM with an hybridization width at the chemical potential much smaller than U. 
Since by the TableHlG' = U, we can safely project out of the low energy subspace the doublet 
(P). The effective AIM within the impurity subspace which includes the singlet (0) and the 
spin-triplet is 

Haim = H^ + Jr {Si + S-^ ■ S 

+ JS,-S2 + ■Jk'^-^T^, (20) 

where and Jk have been defined in Eqs. dH) and (fTIH) . Si and S2 are the impurity 
spin operators for the singly-occupied orbitals 1 and 2, while S and are respectively 
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the conduction-electron spin-density operator and x-component of the pseudo-spin density 
operator, T, at the impurity site. The impurity antiferromagnetic exchange, J = At^/U, 
lowers the energy of the singlet © with respect to the spin triplet. Therefore J alone might 
induce an UFP within the phase diagram, just like in our model as well as in the 2IKM. 
However t± also introduces a scattering potential at the impurity site, last term in the 
right hand side of Eq. ()20|). which is known to be relevant at the UFP. In this respect tj_ plays 
an intriguing role: on one hand it provides a mechanism, the antiferromagnetic exchange 
J, able to stabilize a non-trivial fixed point, but, in the meantime, it also prevents that 
fixed point to be ever reachable. Yet we might wonder whether the critical region around 
the UFP is completely or only partially washed out. In the latter case we should expect 
that the physics of the lattice model close to the MIT is still influenced by the UFP, with 
interesting consequences. We will come back to this issue in Section FVIII 

IV. NUMERICAL RENORMALIZATION GROUP RESULTS 

To study the AIM we used the Wilson Numerical Renormalization Group (NRG) 
method''-^. This technique is known to provide a detailed description of the low energy 
behavior thus allowing a faithful characterization of the fixed points as well as of their 
stability domain. In addition dynamical properties are also accessible by NRG, which, as 
we are going to show, are of notable interest. 

Within NRG the conduction band is logarithmically discretized and mapped onto a one- 
dimensional chain with nearest neighbor hopping integrals which decrease exponentially 
along the chain. The Hamiltonian of a chain with + 1 sites is defined by: 

{N-l 
V A-'^/V (c^ c + c"^ c ) + 

n=0 

+1 (n, - 2) V 2 J [{T^' + {Tyf] I . (21) 

Here U = C^U, J = Ca J, Aq = C| see ©, where Ca = (j^) and all energies are 
measured in units of half the conduction bandwidth. The re-scaling factor A*^^"^)''^ in front 
of (j2T| keeps the lowest energy scale in Hn of order one at each iteration. The original 
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Hamiltonian is recovered in the limit of infinite chain lenght: 

H = lim C;^iA-(^-i)/2^jv (22) 

The size N of the chain determines the infrared cut-off, e.g the temperature T, on a loga- 
rithmic scale (T ~ A^^/^). The method essentially consists in diagonalizing the model on 
a finite size chain, e.g. N, and iteratively increasing the size by one site, from to + 1, 
keeping only the low energy part of the A^-site spectrum. (In our calculations we typically 
kept up to the lowest 2000 states, not counting degeneracies, and used A between 2 and 3. 
We did check that these numbers are sufficient to get accurate results.) 

We restrict our analysis to large values of U, where valence fluctuations on the impurity 
are substantially suppressed. Here, as we discussed, the AIM effectively behaves hke the 
Kondo model (jH)). We fix both U and A and span the phase space by varying the exchange 
parameter J. 

A. Low energy spectrum at the fixed points 

First of all we identify the fixed points by analysing the low energy spectra (with 
typically up to one hundred). Since the conventional size-dependence of the level spacing is 
absorbed by the factor A*^^~^)/^ in front of ()2H) . the low-lying energy levels flow to constant 
values whenever the model is close to a fixed point. Fig. ^ shows that there are two different 
asymptotic regimes separated by a critical value J* < 0. In order to facilitate the interpre- 
tation of that figure, we recall that the ground state of a particle-hole symmetric free-chain 
with + 1 sites is unique if A^ is odd and degenerate if A^ is even. 

For J > J* the low energy spectrum of a chain with odd number, A^ + 1, of sites flows 
towards that of a free chain with an even number of sites and viceversa. This is evident in the 
right panel of Fig. ^where the ground state of the chain with odd A^ becomes asymptotically 
degenerate as for a chain with even A^. Apart from the ground state degeneracy, also the 
low-lying spectrum, i.e. degeneracy and quantum numbers of the levels as well as the level 
spacings, coincides with that of a free chain. As usual, this is as if the first site of the 
chain were locked to form a spin and orbital singlet configuration with the impurity, hence 
becoming unaccessible to the conduction electrons which thus acquire a 7r/2 phase shift per 
conduction channel. It is a conventional Kondo screened phase. 
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FIG. 1: Lowest energy levels versus the chain size N . The right/left panels correspond to a 
deviation 5J / J* = ±4 • 10^"^ from the fixed point value J*. The levels are labeled by the quantum 
numbers {Q,T^ , S), where Q is one-half of the added charge with respect to the ground state value. 

For J < J* the situation is reversed: the low energy spectrum of an odd (even) chain 
flows to that of an odd (even) free chain. Indeed, as shown in the left panel of Fig. ^ the 
ground state with N odd remains non-degenerate for large N. This case corresponds to an 
unscreened phase with the impurity asymptotically decoupled from the conduction bath. 
The phase shift is consequently 6 = 0. 

In between the Kondo screened and unscreened phases we do find a non-trivial fixed 
point, as it is visible in the intermediate cross-over region of the spectrum, see Fig. ^ The 
peculiar non-Fermi liquid character of this intermediate coupling unstable fixed point (UFP) 
is clear by the non-uniform spacing of the low energy levels. A careful analysis of the UFP 
spectrum reveals that it is just the same as that one found in the particle-hole-symmetric 
2IKM4. In Table nil we compare the energies E of the lowest-lying levels of the Wilson chain 
at the UFP, as obtained by NRG, with the prediction x of Conformal Field Theory for the 
2IKM4. The agreement is a clear evidence that the UFP is indeed the same in both models. 

B. Impurity properties at the UFP 

Additional information are provided by the average values of the impurity spin, {S ■ S), 
psudospin, {T -T), and its z-component, ((T^)^). By symmetry, the impurity density matrix 
is diagonal in the six two-electron configurations. The diagonal elements represent the 
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FIG. 2: The angles 6 and as defined through Eq. ((23}. Notice that the fixed point is identified 
by 6* = = 7r/4. 

occupation probabilities P{S, ,T,T^) of states with quantum numbers S, S^, T and T^. 
In the large ?7-limit, where impurity configurations with n ^ 2 have negligible weight, we 
can write 

P(0, 0,1,0) = cos^^, 

P(0, 0, 1, +1) = P(0, 0, 1, -1) = i sin^ 9 sin^ 0, (23) 
P(l, +1,0,0) = P(1,0,0,0) = P(l,-1,0,0) = I sin^^ cos2 0, 

from which it derives that 

(S-S) =2 sin^ e cos2 0, 

(f ■ f) =2 (cos^ e + sin^ e sin^ 0) , (24) 

{{T'f) = sin^^ sin2 0. 

In Fig. 12 we plot the angles 6 and as obtained through (j^^ by the numerical calculated 
average values. The UFP is characterized by ^ = = 7r/4, namely by the value 1/2 of the 
occupation probability of the singlet state The precise value of the other occupation 
probabilities, in other words of 0, are instead not relevant, apart from the obvious fact that 
their sum should be 1/2 too. In fact, if we add the term (fT!?|l with G > 0, we do find the same 
UFP, which locations now depends also on G, which is still identified by P(0, 0, 1, 0) = 1/2, 
i.e. 9 = IT /A, although the weight of the spin-triplet is enhanced with respect to the doublet 
(jni), < vr/4. For large G we do recover the 2IKM values 6 = n/A and = 0, see Fig. 01 
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FIG. 3: The UFP values of 9 and (j) along the path parametrized by the coupling G from our to 
the 2IKM model. 

C. Approach to the fixed points 

As we said the low energy spectrum both in the Kondo screened and unscreened phases 
flows to that of a free chain, with one less site in the former case. The flow towards the 
asymptotic spectrum can be described by a free chain in the presence of a local perturbation 
tern>i^ acting on the first available site, denoted as site 0, of the conduction chain, which 
is actually the second site in the Kondo screened phase. By symmetry considerations this 
local term can be in general written as: 

aa 

+Js* <So ■ <So + Jt* % ■ % — 2 ( + Jj-*) (Xq)^ ■ (25) 

We choose this particular form because it has the advantage that the energy of the center 
of gravity of each multiplet with given charge no is just f/*(r2o — 2)^/2. Upon approaching 
the UFP on both sides, we find that t/* ~ Js* = 7 — > +00, Jt* ~ —67 —00 and 
~ 37/8 ^ +00. The behavior of t^, implies a divergence of the impurity contribution to 
the specific heat coefficient. Namely if 6Cv is the variation of the specific heat with respect 
to its value Cy in the absence of the impurity, then 

-pr~ ~ Pel* 00. 

In reality, it is more convenient to analyse the NRG results by invoking the Fermi liquid 
theory which we present in the Appendix. Through NRG, one can calculate the Wilson 
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ratios related to the conserved quantities'^, namely the total charge, spin and 2;-component 
of T. If 5xi is the variation of the susceptibility with respect to the value Xi without the 
impurity, where Xi = Xc, Xs, Xt charge, spin, and susceptibilities, then the Wilson 

ratios Ri are defined through 

D SXi Cv 

Ki — 



Xi SCv 

On the other hand, Fermi liquid theory implies also that 



Ri 



1 - A;, 



(26) 



where Ai is the dimensionless quasi-particle scattering amplitude in channel i defined in 
Eq. ()A.11|) through the scattering vertex at low incoming and outgoing frequencies and 
the quasi-particle density of states at the chemical potential, see (jA.8j) . In general we can 
introduce a scattering amplitude for each particle-hole and particle-particle channel. In 
particular, besides Ac, As and A^^., we consider the particle-hole scattering amplitudes in 
the channel, which is degenerate with the one, Aj^, as well as in the spin orbital 
channels ST^, ^|t' ST^^'^\ ^st- addition we introduce the amplitudes in the 
particle-particle channels, namely in the spin-triplet orbital-singlet Cooper channel, 
and in the spin-singlet orbital-triplet channels with = and = ±1, respectively. 
As shown in the Appendix, all particle-hole scattering amplitudes can be expressed through 
the particle-particle ones: 

] {6A^ + 2^° + 4Al) , (27) 



Ac 
As 



2A^ - 2^° - AAl) , 
-6A^ - 2^° + AAl) , 
-6A' + 2^0) , 
-2A^ + 2A^o - AAl) , 
-2A^ - 2^0) . 



(28) 
(29) 
(30) 
(31) 
(32) 



Since we are able to calculate by NRG the three Wilson ratios Rc, which is zero in the Kondo 
limit, Rs and R^^, we can also determine the three unknown particle-particle scattering 
amplitudes through Eqs. (pUjl - iP^ . which we plot in Fig. HI 
The first thing to notice is that approaching the UFP, 



•^0 — — ^ST — 



-3, 
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FIG. 4: The scattering amplitudes in the various particle-particle channels as function of J mea- 
sured in units of the UFP J* . Notice the agreement with the values predicted by general arguments 
presented in the Appendix at the UFP, J/\J*\ = —1, at the SU(4) point, J/\J*\ = 0, and in the 
limit of the S = \ two-channel Kondo model, J/\J*\ ^ 1- 

while all the other AiS tend to 1, implying vanishing Wilson ratios. The fixed point seems 
therefore to display a large hidden symmetry, actually an S0(7) as identified by Ref. 3. The 
UFP is equally unstable in the s-wave Cooper channel with S* = 0, T = 1 and = 0, as 
well as in the particle- hole T^*^^^ and ST^ channelaS^i^. All of them correspond to physical 
instabilities as we anticipated and unlike what happens in the 2IKM. On the contrary any 
external field in the other channels do not spoil the UFP, in particular in the charge, spin 
and particle-hole channels, which refer to conserved quantities. 

The physics around and right at the UFP has been uncovered by Conformal Field Theory 
and bosonization^'^'Sii^iik. Due to the existence of two energy scales, the Kondo temperature 
Tk and the exchange splitting, J, the quenching of the impurity degrees of freedom takes 
place in two steps. First, around an energy scale T+ ~ max(TR-, | J|), most of the In 6 
entropy of the two-electron impurity multiplets is removed, leaving behind a residual entropy 
In -\/2 which gets quenched only below a lower energy scale T_ ~ I/7. The latter depends 
quadratically upon the deviation from the UFP, namely T_ ~ | J — J*p/T+. The entropy 
has a low energy linear behavior, S{T) ~ T/T^, followed above T_ by another linear one, 
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S{T) - In V2 ~ T/T+^. At the fixed point, T_ = 0, there is a finite residual entropy 
5(0) = In V2 and S{T) - 5(0) ~ T/T+. A perturbation in any of the relevant channels 
washes out the fixed point cutting off the infrared singularities close to the UFP on an energy 
scale which depends quadratically upon the strength of the perturbation. In Section IVIII 
we analyse more explicitly the stability/instability of the UFP towards symmetry breaking 
fields in particle-hole channels. 



D. Influences of the single-impurity behavior in a DMFT calculation 

Let us now instead discuss the above results in connection with DMFT. Suppose there is 
a lattice model which maps in the limit of large lattice-coordination onto the AIM (0) with 
J < 0. If the model is driven towards a Mott metal-insulator transition, the effective AIM is 
necessarily pushed into a regime in which ~ |</|, namely in the critical region around the 
UFP. As shown in Fig. 01 the s-wave scattering amplitude as well as the equally relevant 
and J^\rp are strongly attractive in an entire interval around the UFP. This suggests 
that the impurity fixed-point instability might transform by DMFT self-consistency into 
a whole pocket where the model generates spontaneously a bulk symmetry-breaking order- 
parameter along one of the relevant channels. As we argued in Ref. [s], if nesting or Van Hove 
singularities are absent, it is most probable that the dominant instability will occur in the 
Cooper channel, the only one which is singular in any dimensions and for any band-structure 
with a finite quasi-particle density of states at the chemical potential. This has been indeed 
confirmed by very recent DMFT calculations in Refs. |^ U- 

The other interesting observation is that in the conventional Hund's regime, the Kondo 
screened phase with J > 0, an attraction in the spin-triplet T = channel develops, < 0. 
In realistic lattice models which map onto the AIM with J > in the limit of large lattice- 
coordination, spin-triplet superconductivity would compete with bulk magnetism. Yet, if 
magnetism is frustrated, spin-triplet superconductivity might emerge. In particular, since 
increasing the Hubbard U in the lattice model implies decreasing Tk in the AIM, which is 
the same as increasing the effective strength of J > 0, we should expect that spin-triplet 
superconductivity is enhanced near the MIT. This has been recently observed by DMFT—. 
However the enhancement of the spin-triplet amplitude is not as dramatic as for the spin- 
singlet one near the UFP at J < 0. This situation would change in the presence of a single-ion 
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anisotropy which favors e.g. spin-triplet pairing with = 0, see the model Eq. (fT^ . As 
we showed, this model is equivalent to upon interchanging the role of T with 5*. This 
suggest that the lattice model which maps by DMFT onto with D > would still enter 
a local critical regime before the MIT. Here it should be dramatically enhanced the tendency 
towards spontaneous generation of a bulk order parameter in the particle-hole channels S^, 
and T S"^ as well as in the spin-triplet Cooper channel with 5^ = 0: ( cL cL — cL 4i ) • 



V. IMPURITY SPECTRAL FUNCTION 



The impurity density of states (DOS), p(e), is defined through 

p(e) = hm [G{e + zv) - G{e - zr^)] , (33) 

where G{ien) is the impurity Green's function in Matsubara frequencies, which, by symmetry, 
is diagonal in spin and orbital indices, and independent upon them. In general 

G{ien)~^ = ien - A(ie„) - S(2e„) = Go{ien)~^ - S(ze„), (34) 

where Go{ien) is the non-interacting, U = J = 0, Green's function, 

k 

is the hybridization function, and E(2e„) the impurity self-energy. Let us suppose to follow 
the behavior of the DOS as the interaction is switched on. We will imagine to increase 
slowly both U and \J\ at fixed U/\J\ ^ 1 with J < 0. When U is small, one can show by 
perturbation theory that 

XmS(e) ~ e^, 

which is the standard result that the quasiparticle decay rate vanishes faster than the fre- 
quency. Therefore at the chemical potential, e = 0, the impurity DOS is not affected by a 
weak interaction, since 

p(0) = lira Xm G {0 + if]) = limlm G o{0 + irj) = — — = pq, (36) 

where Aq = — Xm A(0 + irj) was introduced in Eq. ((7j), and po denotes the non-interacting 
DOS at the chemical potential. In a single-orbital AIM, the above result remains valid 
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even when the interaction is very large. In our case we may expect that something non- 
trivial should instead occur. Indeed, upon increasing U, the AIM enters the Kondo regime, 
with a Kondo temperature exponentially decreasing with U. Therefore at some critical Uc, 
when Tk ~ |</|, the AIM has to cross the non- Fermi liquid UFP. Namely the UFP of our 
AIM can also be attained by increasing the interaction strength, signaling a breakdown of 
the conventional perturbation theory. We now discuss how this criticality shows up in the 
spectral properties. 

The impurity DOS can be obtained by NRG by directly evaluating the spectral function 
Aaaiu;) = i J] |(m|4jn)p5 - (E„ - EJ) (e"^^" + e"^^™) . (37) 

m,n 

For any finite chain A{uj) is a discrete sum of delta-peaks. A smooth DOS is obtained by 
broadening the peaks, which we do following Ref. Q through the transformation 



5{UJ - LOnm) ^ 7 1= ^Xp 



62 



(38) 



where uJnm = En — Em and 6 = 0.55 for A = 2. 

In Fig. El we show the outcome of the numerical calculation. On the Kondo screened side 
of the UFP, the DOS shows a narrow Kondo resonance on top of a broader one. The height at 
the chemical potential is p(0) = po? as expected in a Kondo screened phase. On the contrary, 
in the unscreened side of the UFP, the narrow peak transforms into a narrow pseudo-gap 
within the broad resonance. Numerically we find that p(e) ~ e^. As discussed before, this 
implies that the conventional behavior XmS(e) ~ breaks down across the UFP. Exactly 
at the fixed point, both the narrow peak and the pseudo-gap disappear, leaving aside only 
the broad resonance. The calculated DOS at the chemical potential seems to be half of 
its non-interacting value, see Fig. |S| In other words our numerical results point to a DOS 
at the chemical potential which jumps across the UFP, being p(0) = po everywhere in the 
Kondo-screened phase, p(0) = in the unscreened one, and p(0) = po/2 right at the UFP. 



VI. MODELING THE IMPURITY DENSITY OF STATES 

It is possible to infer an analytical expression of the impurity DOS. First of all we notice 
that the values at the chemical potential in the screened and in the unscreened Kondo 
regimes are compatible with general scattering theory. In both phases the impurity has 
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FIG. 5: Impurity DOS in the presence of particle-hole symmetry across the fixed point. The 
temperature is set by the lenght of the chain; it is practically zero. In the upper panel we draw the 
DOS's well inside the Kondo screened phase {J/ J* = 0) and the unscreened one {J /J* = 5.75). 
Here U = 2, Aq = U /(67r) and J* turns out to be ~ —0.0035, all in units of half the conduction 
bandwidth. Notice the narrow peak which transforms into a narrow pseudo-gap. In the lower panel 
we show in more detail the behavior of the low energy DOS across the UFP. (From top to bottom 
J/ J* = 0, 0.859, 0.945, 0.988, 1.002, 1.031, 1.146) 



disappeared at low energy, either because it has been absorbed by the conduction sea or 
because J has taken care of quenching the impurity spin and orbital degrees of freedom. 
This in turns means that what remains at low energy is just a potential scattering felt by the 
conduction electrons plus a local electron-electron interaction term. The on-shell 5'-matrix 
at the chemical potential has in general elastic and inelastic contributions (see Ref. llTI). At 
zero temperature only the former survives. Since we considered just s-wave scattering, the 
elastic component of the S'-matrix is given by: 



5(0) = 1 - 27iipcT{0) = 1 - 2 TT Ao p(0), 



(39) 
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where pc is the conduction electron DOS at the chemical potential per spin and band, and 
the T-matrix is defined through the conduction electron Green's function Q by 



Q = Qo + GqT Go- 
On the other hand the S'-matrix is related to the scattering phase shift by 

^(0) = e2^^W. (40) 

In the Kondo screened phase, we know that 6{0) = 7t/2 which, through pUj) and ()39|) implies 
p(0) = I/ttAq, namely its non-interacting value po- On the other hand, in the unscreened 
regime 5(0) = hence p(0) = 0, as we indeed find. It has been proposed that at the non- 
Fermi liquid fixed point of the overscreened S = 1/2 two-channel Kondo model the S'-matrix 
is instead purely inelastio^'^SiiS. That would imply a vanishing elastic contribution, 5'(0) = 
in dnni), and in turns a DOS at the UFP 

P(0) = 7^ = Ipo, (41) 

which is indeed compatible with our numerical resulta^S. Yet there is a difference between 
the UFP of our model, equivalently of the 2IKM, and the non-Fermi liquid fixed point of 
the S = 1/2 two-channel Kondo model. While in the latter the specific heat has a singular 
temperature behavior right at the UFP, in our model it has a conventional linear behavior. 
The above observation suggests the following simple analytical expression of the low-energy 
impurity DOS: 

^*W-f(^±J^). (42) 

where the plus sign refers to the Kondo screened phase and the minus to the unscreened 
one. The two energy scales have the same meaning as in the previous Section. In particular 
T_ controls the deviations from the UFP, so that right at the UFP, when T_ = 0, the DOS 
is 

The model-DOS (jl^ also implies a model impurity Green's function in Matsubara frequen- 
cies: 

G±(«e„) = ^ ± ^ 1 . (44) 

2 Ao V + u + signe„ «e„ + «T_ signe„ / 
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FIG. 6: Fit values of and T_ in the Kondo screened phase. The hne is a quadratic fit, 
T_ = A{6J)^. The Hamiltonian parameters have the same values as in Fig. El 

The fixed point Green's function, G*(ie„), is identified by T_ = 0. The impurity self-energy 
can then be extracted by the relation 



In particular, at low frequency we find that 



An / 1 1 



(45) 



in the Kondo screened phase, hence a standard linear behavior. On the contrary, in the 
unscreened regime the self-energy is singular 

1 2AoT+T_ 



ze„-S_(ze„)- . 

ten T+ - T_ 

Finally, at the fixed point the self-energy is finite at zero frequency, being given by 



(46) 



ien - SJze„) = iAo 



T. 



(47) 



+ 



We have checked that the model-self-energy gives indeed a good representation of the actual 
numerical results. In Fig. El we draw the fit values of T+ and T_ within the Kondo screened 
phase. 

We can further test the consistency of the approach by invoking the scattering theory 
which, by the Friedel's sum rule, allows us to identify the scattering phase shifts through: 



5(e) = ImlnG'(e + «0^ 



(48) 
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By means of our ansatz for the impurity Green's function ()44|) we readily find that the 
expression of the low-energy phase-shifts is 

+ (49) 

within the Kondo screened regime, and 

(^ + ^) ^«_6, (50) 

in the pseudo-gap unscreened phase, consistent with our starting assumption. Moreover, by 
the energy dependence of the phase shifts, we can calculate the impurity correction to the 
specific heat 

VII. PARTICLE-HOLE SYMMETRY BREAKING TERMS 

In this Section we analyse more in detail various symmetry breaking terms in the particle- 
hole channel. In particular we are going to consider the three following perturbations to the 
original Hamiltonian (Q) with z/ = 0: 

5Hp^h = vUnd = Ud, (52) 
= h,T\ (53) 
= KT^. (54) 

The term ()52|) breaks particle-hole symmetry trying to occupy the impurity with 2 — u 
electrons instead of two, see (PJ). The other terms, and (jH^ . split the orbital degeneracy. 
Among them, only 6Hx is predicted to be relevant and wash out the fixed point, at least 
according to bosonization^. Actually this looks a bit strange result if one invokes naively the 
argument of Ref.y to demonstrate the existence of an UFP in the absence of any particle- hole 
symmetry breaking term. This argument is based on the observation that, when 0{2)orbit 
symmetry holds, the phase shifts in both orbital channels have to be equal, 6i = 62- By 
general particle-hole symmetry, this further implies that 26i = 262 = mod{Tx). Since for 
J ^ Tk > we know that 5i = 82 = 7r/2, while for J <^ —Tk < 0, 5i = ^2 = 0, there must 
necessarily be a fixed point in between. 
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Let us assume now that the T^-term fl53|) is present and follow Ref. y to demonstrate 
that the necessary condition for the existence of an intermediate fixed point does not hold 
anymore. Since (jK^ is absent, there is still a residual particle-hole symmetry according to 
which 

^1 + ^2 = mod{T[). 

If 5i = —62 then the two limiting cases, 61 = 62 = and 61 = —62 = 7r/2, can be smoothly 
connected without requiring any critical point in between. This argument thus proves that 
an intermediate fixed point does not need to exist, yet it does not demonstrate its non- 
existence. Indeed we know by bosonization and we now show by NRG that both as 
well as ()52j) do not wash out the UFP. On the contrary a T^-term fl54|) does destabilize the 
fixed point, as shown later. 

Let us go back to Eq. (jH^ and try to guess how would it change in the presence of flS^ 
and/or We have now to introduce one S'-matrix for each channel, Sa with a = 1,2, 

satisfying 

-JZeSaiO) = cos 26a{0) = 1 - 2 7r AoPa(O). (55) 

Let us assume that, across the UFP, the zero-frequency phase shifts still jump by 7r/2. In 
other words, if we denote as 

5-,a(0) = 6a, (56) 

the phase shift in the unscreened phase, in the Kondo-screened one the phase shift should 
be 

5+,„(0) = 5a + |. 

Through ()55p this would imply a jump of the DOS at the chemical potential given by 

p+,a(0) - p~,a{0) = COS 25a = Po cos26a. (57) 

The above scenario predicts that although the pseudo-gap in the unscreened phase is partly 
filled away from particle-hole symmetry, yet the DOS has a finite jump across the UFP. This 
is indeed confirmed by NRG. In Fig. Owe plot the DOS at fixed u = 0.05, see Eq. (jK^ . 
across the UFP, clearly showing the jump. We notice that if only ()52|1 is present, then 
6i = 62 in Eq. If (jK^ is absent but (jK!^ is present, then 61 = —62, yet the behavior 

across the UFP is similar, which is the reason why we just show the results with finite u. 
This behavior is also compatible with the NRG result that the charge and Wilson ratios 
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FIG. 7: Impurity DOS across the UFP in the presence of a finite v = 0.05 which breaks particle- 
hole symmetry. From top to bottom J/ J* = 0,0.28,0.57,0.86,1.14,1.43,1.71. Notice that the 
DOS at the chemical potential is always finite, although very small hence not visible in the figure. 
U and Aq have the same values as in Fig. 13 

vanish around the UFP. Actually they all suggest that the model can absorb a chemical 
potential shift, equal or different in the two channels 1 and 2, on a high energy scale, at 
least of order T+, without having to modify what takes place at lower energies of order T_: 
a kind of Anderson's compensation principle for our conserved quantities. Following these 
observations, we argue that the DOS for orbital a = 1,2 in the presence of any of the two 
perturbations, and (|^. assumed to be weak, can be modeled as 



P±,a(e) = — 



± COS 2(5(j 



(58) 



where again the plus refers to the Kondo screened phase, the minus to the unscreened one. 
Pa = p+,a(0) is the value of the DOS at the chemical potential in the screened regime, while 

f^±,a = ±T+ sin2(5a. 

According to the model DOS (j3Hj) . the narrow peak and pseudo-gap remain pinned at the 
chemical potential, e = 0, while only the broad resonance moves away from particle-hole 
symmetry. 

Let us now study what happens if, starting from the particle-hole symmetric pseudo-gap 
phase we move away by increasing z/, keeping all other Hamiltonian parameters fixed. As 
shown in Fig. |H1 z/ is able to drive the model across the UFP. This result could be foreseen. 
Indeed u forces the impurity to accomodate 2 — z/ electrons. If u = 1, the impurity tends to 
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FIG. 8: Impurity DOS upon increasing the strength of particle-hole symmetry breaking v starting 
from the unscreened pseudo-gapped phase {y = 0,0.05,0.1) up to the Kondo screened one (i^ = 
0.15,0.2). In the left inset it is shown the low energy part across the UFP (from top to bottom 
= 0.2,0.175,0.15,0.125,0.1,0.075,0.05,0); notice the analogy with the p-h symmetric case in 
Fig. [SJ In the right inset we explicitly show the gradual filling of the pseudo-gap upon increasing 
z^. The values of U and Aq are those of Fig. [3 

be singly-occupied. Therefore in the Kondo limit it behaves like a spin S = 1/2 and pseudo- 
spin T = 1/2 moment, which can be perfectly Kondo-screened and it is moreover stable with 
respect to little changes of u with respect to = 1. Hence, if the model is at = in the 
pseudo-gap phase, it has to cross a fixed point to reach the large-z/ Kondo screened regime. 
This behavior is quite interesting in connection with DMFT lattice calculations, since it 
implies that the lattice-model local critical regime, which reflects the single-impurity UFP, 
may also be attained by doping, as recently confirmed^. 

A completely different behavior occurs if we introduce instead a perturbation of the 
form ()54|1 . Here, as expected, we do not find any jump of the DOS, as clear in Fig. 121 where 
we compare the DOS at the chemical potential in the presence either of (jK^ . hp^h 7^ 0, 
or (jnH), hx 7^ 0. This demonstrates that a perturbation in the particle-hole channel which 
breaks the orbital 0(2) symmetry is relevant at the UFP, unlike fl52p and which instead 
do not break the 0{2) orbu symmetry. 

Finally let us discuss what happens in the AIM which corresponds within DMFT to two 
Hubbard planes coupled by a transverse hopping, Eq. ()18|1 with J = V = plus the term 
(I19|) . We already noticed that t± plays an ambiguous role. It generates an antiferromagnetic 
exchange, J = At^/U, which may stabilize an UFP, but it also induces a relevant 
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FIG. 9: Comparison of the DOS values at the chemical potential as function of J either in the 
presence of a finite particle-hole symmetry breaking hp-h, p{0, hp^h), or of a symmetry breaking 
hx, p{0, hx), normalized to their values at J = 0. Notice that hx, although three order of magnitude 
smaller than /ip_h, washes out the DOS jump contrary to hp-^h- 

perturbation, see Eq. ((201). Since the UFP is never reachable, the model always flows to a 
Fermi liquid fixed point. In the presence of it is more appropriate to introduce the even 
and odd combinations of the orbitals 1 and 2: 

dea = {di„ + d2a) , 

and correspondingly the even and odd conduction-electron scattering channels. According 
to what we said at the beginning of this Section, we expect the phase shifts (5e = —bo to be 
smooth functions of J. If there were no remnant of the UFP, the DOS's should simply show 
a resonance, the even channel above the chemical potential and the odd channel below it. In 
reality the behavior of the DOS remains strongly influenced by the UFP, even though never 
reachable. This is evident in Fig. [TUl where we draw the DOS of dgo-, Peis)-, (the odd one is 
simply obtained by reflection around zero energy) at fixed t_L upon varying the hybridization 
width Aq. There is no point at which the DOS jumps at the chemical potential, yet a partly 
filled asymmetric pseudo-gap remains. In Fig.^Jwe draw the low-energy difference between 
the even and odd DOS's, which is also the off-diagonal spectral function Ayiis). Ayi^e) shows 
a low energy feature which has a non-monotonic behavior in Aq and almost develops into a 
singularity around Aq = 0.47. We think that these results bring to the fore that alone is 
able to drive the model very close to the UFP. In other words the width of the critical region 
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FIG. 10: Impurity DOS of dea, Pe{^), for the AIM Eq. H2U() . The different curves correspond from 
the top to the bottom to values of Aq = 0.5, 0.47, 0.45, 0.4, 0.3 with t± = 0.05 and U = 8. These 
values correspond to Jk = 0.08, 0.075, 0.072, 0.064, 0.049 and J = 4:t'j_/U = 0.00125. We notice the 
remnant of an asymmetric pseudo-gap of order J. 

is larger than the energy scale which cut-off the fixed-point singularities, although both are 
generated by the same t±. 



VIII. CONCLUSIONS 



In this work we have analysed the spectral properties of the two-orbital Anderson impurity 
model, Eq. which includes an exchange splitting J which favors, if negative, a non- 
degenerate impurity configuration. This model was already shown in Ref. '.8 to posses a 
non-Fermi liquid fixed point which separates a phase where conventional Kondo screening 
takes place from an unscreened phase in which J takes care of quenching the impurity 
degrees of freedom. 

The impurity density of states has the following behavior across the fixed point in the 
presence of particle-hole symmetry. In the Kondo screened phase it displays a conventional 
very narrow Kondo resonance on top of a broader resonance. On the contrary, in the 
unscreened phase a narrow pseudo-gap appears within the broad resonance. At the fixed 
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FIG. 11: Off-diagonal spectral function, Ai2{e), with t± = 0.05. The different solid curves corre- 
spond from the top to the bottom for e > to values of Aq = 0.47, 0.45, 0.4, 0.3, while the dashed 
curve corresponds to Aq = 0.5. We notice that the low energy feature first moves towards zero 
energy when Aq increases from 0.3 to 0.47, but from 0.47 to 0.5 it goes back again. Moreover, 
around Aq = 0.47, vli2(e) is almost singular. 

point only the latter survives. Away from half-filling, the pseudo-gap remains pinned at the 
chemical potential, although it gets partly filled. Yet there is still a fixed point across which 
the density of states at the chemical potential jumps. Finally we have explicitly shown that 
the intermediate fixed point is unstable towards physical symmetry breaking fields, which 
include both particle-hole and particle-particle channels. The relevance of this impurity 
model for Dynamical Mean Field theory calculations has been already emphasized in Ref. la 



and confirmed by Refs. IflUlll. Here we would like to clarify some aspects in view of the newly 



discovered spectral properties. 

As discussed in Ref. ^ any lattice model which maps by DMFT onto the impurity model 
plus ()13|) should encounter the unstable fixed point before the Mott transition, namely 
when the effective quasiparticle bandwidth becomes of the order of \J\. However the in- 
stability of the single-impurity fixed point should likely transform into a bulk instability 
through the DMFT self-consistency conditions. As we showed there are several compet- 
ing physical instabilities around the fixed point, in the particle-hole and particle-particle 



31 



channels. In the absence of nesting or van Hove singularities, we argued in Ref. |8| that 
the particle-particle channel dominates, leading to a superconducting pocket just before 
the Mott transition, which has been indeed observed by DMFT^. However there might be 
physically relevant cases where those band-structure singularities occur, which would fa- 
vor a uniform or modulated order-parameter in one of the particle-hole unstable channels. 
What should we expect upon moving away from these peculiar cases, for instance by dop- 
ing? Clearly the band-structure singularities weakens upon doping. Yet the fixed point is 
not washed out away from particle-hole symmetry. We showed in fact that the pseudo-gap 
remains pinned at the chemical potential. We believe that this would result into a compe- 
tition between particle-hole and particle-particle channels which gradually turns in favor of 
the latter, thus predicting a particle-hole order parameter which dies out upon doping in 
favor of a superconducting one. 

Equally interesting is what we find for the Anderson impurity model which corresponds 
within DMFT to two Hubbard planes, with large in-plane coordination, coupled by an 
hopping term t±, Eq. ()18|1 with J = V = plus the term ^T^i . Here the physics is not 
as transparent as in the model (0), essentially because t± provides at the same time a 
mechanism for the existence of a fixed point as well as for its instability. However the 
numerical renormalization group results for the single- impurity show evidence that an almost 
critical region does exist, in spite of the fact that the non-trivial fixed point can never be 
attained. This suggests that the physics of the two coupled Hubbard planes close to the 
Mott transition may still be influenced by the single-impurity fixed point. 

Finally, we briefly comment what our results would imply for the model with conventional 
Hund's rules, see Eq. (fT^ . This case in the absence of single- ion anisotropy corresponds 
to a Kondo screened regime where, as we showed in Fig. HI the spin-triplet Cooper channel 
is attractive. By increasing the Hubbard U in the lattice model, the Kondo temperature 
of the effective AIM decreases, which implies that the spin-triplet dimensionless scattering 
amplitude gradually also decreases, see Fig. IH for increasing J > 0. This suggests that 
the instability towards spin-triplet superconductivity may actually be enhanced by strong 
correlations, compatible with recent DMFT calculations^^. In addition we would expect 
that, in the presence of a single- ion anisotropy, D > in Eq. (fTBj) . the enhancement of 
spin-triplet superconductivity might be even more dramatic. 
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APPENDIX: FERMI LIQUID THEORY OF THE ANDERSON IMPURITY 
MODEL 



In this Appendix we build up a Fermi liquid theory of our AIM closely following the 
conventional approach (see for instance Ref. l21l). Our purpose is twofold. First the Fermi 
liquid theory provides a framework to analyse the NRG data. Moreover it allows to introduce 
within DMFT the concept of a local Fermi liquid description in addition to the conventional 
one, which refers instead to low frequency and momentum scattering amplitudes. 

Let us consider more generally a multi-orbital Anderson impurity model. As in our case, 
we assume that besides spin rotational symmetry also orbital degeneracy is preserved, so 
that the fully interacting impurity Green's functions are diagonal and independent either 
upon spin and orbital indices. 

The variation of the electron number with orbital symmetry a and spin a associated with 
the presence of the impurity is given by^ 



where the integration contour encloses clockwise the real axis, f{z) is the Fermi distribution 
function in the complex plane and Gaa the impurity single-particle Green's function. Since 
the Green's function has a branch cut on the real axis, the above expression is also equal to 



with 6 an infinitesimal positive number. The impurity density of states is further determined 
through 



An, 



aa 





(A.l) 




(A.2) 



If we introduce a source field in the Hamiltonian by 




aa 



where 




k 



then 




1 



'^aa 
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where 

Aaaiiu^n, Ka) = iHc^ ^ , , , 

k 

is the hybridization function in the presence of the source. Therefore the derivative with 
respect to the external field of the variation of the electron number associated with the 
impurity is given by 



K=e+isJ V /fe=0 



where 'Eaa{i^n) is the impurity self-energy and we made use of 

f dAaa{z,ha^) \ _^ . dAjz) 
\ ~ 'Jab'Ja/S — ^ , 

V dhbfs dz 
being A{z) the hybridization function in the absence of h. On the other hand 

V oiT'bfi Jh=o f^^^ V ofze^ J 

(A.4) 

where we used the property that, at /i = 0, the Green's function does not depend on a and 
a. The interaction vertex is the reducible one. 

Let us assume that there exists a set of conserved operators 

k abafi 

+Y.dL{m:Ub,^ 

abap 

where M*^*-* are hermitean matrices and the suffix i identifies the particular conserved oper- 
ator. For convenience we adopt the normalization Tr ^M*^*^ • M*^*^^ = 1. Then, if we add a 
source field 

we can use the basis which diagonalizes M*^*) and apply the above results to find the variation 
of (tM '-*-') associated with the presence of the impurity at first order in the applied field. 
Going back to the original basis, we would find the following expression of the difference 



35 



(5x^'^ between the susceptibilities in the presence and absence of the impurity: 

^ [ dh») , 



h=0 



de df{e) 
TT de 



Im I G{e + iS) 



1 - 



dA{ie] 
die 



ie=e+iS 
■yS 



7^ X] X] r6^,d(5;c7,oa(e + i^, i^n] ^^n, ^ + i^) 

n abed a0^5 

dA{ierr, 



dier 



(A.5) 



Hereafter we drop the suffix i. One can demonstrate that the following Ward identities 
hold for the impurity 

[E(i6 + ia;)-E(i6)]M„f 



~^ X XI raa,d5;c7,b/3(^e + ^ i^n, i^n + , «e) 
n cd;75 



(A.6) 



It then follows that 



^^^Kb = -^XXX^«°''^'5;^'^'3(^^'^^"'^^"'^^)^^ ^(^^")' 



- lim ^ r„c,,(i5;c7,6/3(«e + iuj, ien] ien + iuj, ie) 



M]^G{ien + iu;)G{ien) 



[-A{ien + iuj) + A(ien) 



dA{ie„ 
dien 



--^ X X raa,d5;c7,&/3(^e, ^e„; ie„, ie) M^f G{ienf i 1 

n cd ■yS ^ 

/ r„„,d5;c7,6/3(«e, 6 - z5 ; e + zf!) , ze) 

zn ae ^^.^^ 



(A.7) 



Let us define the quantity 



p-. = |^-^Im|G(e + ^^) 



1 - 



9A(ie) 
die 



dE{ie) 
die 



(A.8) 



36 



which plays the role of the quasiparticle DOS at the chemical potential. Then, through 
lU, ()A.7|) and ()A.8|) . the following equation is readily found 



= EE^. (*'")r 

ab a/B 



^aa,dS;n,bi3{(^ + iS, e' - iS'; e' + i6',e + i6) 

cd;7(5 



, , , dm am 

cd;7(5 



(A.9) 



The last expression is obtained by noticing that only the imaginary part of G{e + i6) con- 
tributes, where Im G{e+i6) = —npi^e). Eq. ()A.9j) allows to express any susceptibihty to fields 
coupled to conserved quantities. If the hybridization function is smooth at low energies, then 

A{e' - i5') - A{e' + iS') - 2iAo, 



hence we can rewrite (IA.9I) as follows 



Ao_ 

p*7r 



de 



dt' 



= [1 - Ai] , 

which allows to identify local Landau A-parameters through 



(A.IO) 



A, 



de 



de' 



rac,,d5;c7,b/3(e + ^^5, c' " ^5'; c' + iS' , e + 15) 
{M^^r:d{M^')b:G{e' + ^6')Gie'-^6')]. 



(A.ll) 
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The above expression is quite general but simplifies substantially when the imaginary part 
of the impurity self-energy vanishes at low real frequency. In this case 

1 



where = e^^"* + ReS(O) is the actual position of the (i- resonance. Then, through ()A.2I 

1 Aq a 

.2 I A 2 

Analogously 



p(0) = -3XX2 = -^G{^Q^)G{^Q-)■ (A.12) 
+ A^ TT 



_ p(0) 1 _ ^ fdn^e) 



Pa / ^ ,.,\'yS 



Z ' Z \ die 

hence 

abed 0/375 

which is the more conventional expression of the Landau parameters^!. Although the above 
equation is a particular case of the general one ()A.11|) . to simplify the notations in what 
follows we will use ()A.13|) as a short-hand expression of ()A.11|) . 



1. Application to the twofold orbitally degenerate AIM 

Let us now apply the above results to our model. An incoming pair can be a spin-triplet 
orbital-singlet, with a scattering vertex at zero incoming and outgoing frequencies given by 

ri(T,2o-;2(T,lo-, -^la,2-i7;2 -17,1(7 "~ la,2 -a;! -(7,2a ■ 

Here 1 and 2 label the two orbitals with = +1/2 and = —1/2, respectively. Alterna- 
tively it can be a spin-singlet orbital-triplet with = 0, with scattering vertex 

-i- 1(7,2 -a;2-a,la "T -- 

or with = ±1, in which case 



Tq ~^ 7:^l(7,2-a;2-a,la + la,2 -a;l -a,2(T , 



r± — ^ Tier,! -(t;1 -o-.lcr, ^ 2a,2 -(7;2 -(7,2a ■ 

In reality it is more convenient to introduce the dimensionless scattering vertices: 

= Z^p.Ti = Zp{0)T\ 

At = Z'p.r'o = Zp{0)T'o, (A-14) 
Al = Z^p,Tl = Zp{0)Tl. 
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As we previously showed, only the susceptibilities of conserved quantities can be expressed 
in terms of the Landau parameters ()A.13|) . which are simply connected with the scattering 
vertices at zero frequency. Yet we can still define Landau parameters for non-conserved 
quantities, which, although do not serve to calculate susceptibilities, may provide a qualita- 
tive estimate of their magnitude. Therefore we are going to introduce the Landau parameters 
for the charge, Ac, spin As, the 2;-component of the pseudo-spin T, A^^, all being related 
to conserved quantities, but also for the x and y components of T, Aj^, as well as for the 
spin-orbital components, A^gj, and A^j,. In terms of the dimensionless amplitudes ()A.14|) 
they can be shown, after some lengthy algebra, to have the following expressions: 



Ac 
As 

y4.!j-t 



A 



ST 



A 



ST 



{2A^ - 2^° - 4^i) , 
(-6^^ - + 4^0 ) 

{-2A^ + 2^° - 4^^) 



(A.15) 
(A.16) 
(A.17) 
(A.18) 
(A.19) 
(A.20) 



Let us consider several possible cases. 

• If J = 0, SU(4) symmetry holds. Then = A^q = A^^ = A, leading to 



Ac 
As 



3 A 



-A. 



In the s-d limit, when the AIM maps onto an SU(4) Kondo model, the charge com- 
pressibility in negligible, leading to 3 A = 1. The Wilson ratios for the conserved 
quantities are defined trough 

5x'^^ Cv 



Ri 



l-Ai, 



Xo 5Cv 

where ^x''*'' has been defined in ()A.5|) . Xo = Pc and Cy are respectively the conduction- 
electron susceptibihty and specific heat in the absence of the impurity, and 

xr^ P* 
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is the variation of the specific heat due to the impurity. Hence all Wilson ratios have 
an universal value, 

Rs = R^ = RsT = l- A = 4/3, (A.21) 
in agreement with Conformal Field Theory. 

If J ^ Tk > the impurity gets frozen in the Kondo limit into a spin S=l. Then 
both Ac = 1 and A^^ = 1, which implies 

Al = 1, 

However one expects that, being the spin-triplet an orbital singlet, the SU(2) orbital 
symmetry gets restored at the fixed point, much in the same way as spin anisotropy 
is irrelevant at the Kondo fixed point. This further implies that 

namely As = —5/3, with a Wilson ratio Rs = 8/3, in agreement with known results. 

Let us now suppose to be close to the UFP within the Kondo screened regime. As 
usual the charge degrees of freedom are suppressed already below U, so that we can 
still assume Ac = 1. Moreover we expect that the spin and the orbital degrees of 
freedom related to get quenched below T+, while the remaining ones only below 
T_ ^ T+. Therefore at very low temperatures T < T_, we can safely assume that 



T_ 5xs ~ T_ 5xt ~ 77^ 0' 



namely As = Aj^ = 1. As a result we find that 

Al = A^ = 1, (A.22) 
A^o = -3. (A.23) 

Eq. (|X:23|l impl ies a strongly attractive s-wave singlet channel. The other Landau 
parameters are thus given by 

A^ = aIt = -3, (A.24) 

Ajj. = 1. (A.25) 
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This further proves that the fixed point is equally unstable in the s-wave Cooper 
channel Fq, as well as in the T^, and ST'^ particle-hole channels. 

We finally notice that, although the Landau A-parameters would suggest that the 
susceptibilities in the unstable channels, all of which correspond to non-conserved 
quantities, diverge as 1/T_, in reality they only diverge logarithmically^ii^^. This is 
not incompatible with Fermi liquid theory, which allows to express in terms of the 
A-parameters only those response functions related to conserved quantities. 

Let us now use our model self-energy to extract some additional information. Through 
Eq. we find that in the Kondo screened regime the expression ()A.13|) holds with a 

quasiparticle residue 

Indeed Z ~ 2r_/Ao upon approaching the unstable fixed point. 

On the contrary, the general expression ()A.11|) has to be used inside the non-Kondo 
screened pseudo-gap phase. Through Eq. (ji^ for we find that at low frequency 

^ ^ ^ 4A2 TlT^_ 2Ao ^ ^ r+ + T_ 

By Eq. ipUl) the quasiparticle DOS at the chemical potential turns out to be finite, 

p. = -^1^, A.27 

even though the impurity DOS vanishes. In conclusion, within the pseudo-gap phase the 
Landau parameters have the following expression 

2 (T+ + T_)' Loo de ^ dt' ^ 

(M^y)i (A/<'))*] . (A.28) 

In spite of the anomalous impurity Green's function, the low-energy behavior should still be 
described within a local Fermi liquid scenario by finite Landau parameters Aj's. Therefore, 
since the impurity DOS vanishes quadratically in the pseudo-gap phase, then the scattering 
vertices must display a singular behavior 

rCe e'- e' e] ^- 

^ ' ' ' ^ (e + e')^' 
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to compensate for the vanishing DOS's and provide finite A's. 
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TABLE II: Energies E of the low energy levels and their degeneracy deg at the unstable fixed point. 
The levels are labeled by the quantum numbers Q, half of the deviation of the number of electrons 
with respect to the ground state, S, total spin, and T^, total z-component of the pseudo-spin. The 
value X is the prediction of Conformal Field Theory for the two-impurity Kondo model^. Notice 
the anomaly of the member within the (1/2, 1/2, 1/2) multiplets identified by a *, which was also 
found in Ref. J. There an explanation for the discrepancy was proposed. 
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